Microscopic Model for Granular Stratification and Segregation 
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We study segregation and stratification of mixtures of grains differing in size, shape and material 
properties poured in two-dimensional silos using a microscopic lattice model for surface flows of 
grains. The model incorporates the dissipation of energy in collisions between rolling and static 
grains and an energy barrier describing the geometrical asperities of the grains. We study the phase 
diagram of the different morphologies predicted by the model as a function of the two parameters. 
We find regions of segregation and stratification, in agreement with experimental finding, as well as 
a region of total mixing. 

PACS Numbers: 05.40+j, 46.10+z, 64.75+g 

When mixtures of grains of different sizes are 

poured on a heap, a size segregation of the mixture is 
observed; the large grains are more likely to be found 
near the base, while the small grains are more likely to 
be near the top [p|-|T^]. If the grains differ not only in 
size but also in shape and roughness, a spontaneous pe- 
riodic pattern arises upon pouring the mixture in a two- 
dimensional cell. When a mixture of large-cubic grains 
and small-rounded grains is poured in a vertical Hele- 
Shaw cell (two vertical slabs separated by a gap of ap- 
proximately 5 mm) the mixture spontaneously stratifies 
into alternating layers of small-rounded and large-cubic 
grains Jl5| . Otherwise, the mixture only segregates when 
the large grains are more rounded than the small grains 
|l5f with the large-rounded grains being found near the 
bottom of the cell. 

The dynamical process leading to stratification was re- 
cently studied numerically and theoretically fTfj| ], using 
the set of continuum equations for surface flows of gran- 
ular mixtures developed in 17-23]. The physical quan- 
tities defined in this phenomenological formalism are to 
be understood as an average over a certain coarse grain 
length on the surface of the sandpile (larger than the size 
of the grains), where hydrodynamic equations are valid, 
and any quantity defined below this scale is not well- 
defined. Thus the relevant length scale appearing in this 
formalism is that of the coarse grain scale (ss 5d) and 
not of the grain size (d). However fluctuations may also 
occurs at the level of the grains, so that a microscopic 
description of collisions and transport of grains may be 
needed to describe this situation. In this paper, we study 
the dynamical segregation process in two-dimensional si- 
los by using a microscopic model of grain interactions. 

We start by defining the model for the case of a single 
species pile pl| , p2[ , and then we discuss the generaliza- 
tion to two types of grains differing in size and shape. 
The microscopic model is defined on a square lattice. 
Each grain has width and height of one pixel. Since the 
experiments are done by pouring a fixed flux of grains, 
we deposit N grains at a given time step at the top of the 



first column of the pile. The grains start with a certain 
initial kinetic energy eo, which will be lost in collisions 
with the static grains of the pile as the grains move down 
the slope. Only one rolling grain is in contact with the 
pile surface, and therefore interacts with the static grains 
of the pile. The remaining rolling grains are convected 
downward with unit velocity without loosing their energy, 
i.e., they move to the nearest-neighbor right column. The 
loss of energy of the rolling grain interacting with the 
surface pile is determined by the restitution coefficient, 
r, which gives the loss of energy per unit time. The in- 
teracting grain moves until its energy is smaller than a 
certain energy barrier u and stops. When the interacting 
rolling grain stops, one of the remaining rolling grains 
starts to interact with the surface pile until it losses its 
energy and stops. When all the N rolling grains stop, a 
new set of N grains is dropped at the first column of the 
pile, and the same rules are applied again. 

The dynamics of a rolling grain with energy e interact- 
ing with a static grain located at height h% are defined as 
follows: 

• We test if the nearest neighbor position i+1 is ener- 
getically favorable, by first calculating the energy 
test etest, defined as the energy the grain would 
have after moving to the new position i + 1, ac- 
cording to 



where Ah 



etest = (e + Ah) * r, 



hi 



(1) 



Then, the energy test is compared with the energy 
barrier u. If etest > u, the grain moves to the near- 
est neighbor position i + 1, and the energy of the 
grains is updated: e = etest- Then the procedure 
is applied to the interacting grain at column i + 1, 
and so on. If etest < u, the interacting grain stops 
at position i and increases the height of the pile at 
i by unity. The remaining rolling grains are always 
convected downward to the position i+1. 
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Simulations and analytical calculations |21| show that 
the angle of repose of the pile is an increasing function 
of the energy barrier u, and a decreasing function of the 
restitution coefficient r. 



721 < "u < #22 < 



(2) 
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FIG. 1. Phase diagram for surface flows of granular mix- 
tures predicted by the present model. In all simulations 
we deposit two grains of type 1 and two grains of type 2 
(TVi — N2 = 2), and we set eo = 1, and i[> u — 0.3. For the 
region where S r > we use S u — 0.3, for the region where 
6 r < —0.3, and for the intermediate region of 

S r — we use S u — 0. The letters indicate the simulations 
shown in Figs. |i| and []. 

Next we generalize the model to the case of two types 
of grains. In stratification was reproduced using a 
discrete model defining different critical angles between 
the grains. The critical angle is defined as the maximum 
angle at which a rolling grain will be converted into static 
grain. This angle depends on the type of rolling grain and 
the type of static grain which is interacting with. Thus, 
for two types of grains there are four different critical 
angles also called generalized angles of repose 9 a p, with 
a,/3 = 1,2. Stratification of grains differing in size and 
shape is the result of a competition between size seg- 
regation and shape segregation fl23fl. This competition 
was incorporated in the models of ]16| at a macroscopic 
level, by considering certain relations between the gener- 
alized angles of repose. The angle of repose of the pure 
species depends on the shape of the grains: the rougher 
the grains the larger the angle of repose. Thus, for mix- 
tures of cubic grains (type 2) and rounded grains (type 1) 
we have #22 > #11 • On the other hand, if the grains have 
different size, the cross- angles of repose 9 a p are different. 
Since small grains roll down on top of large grains easier 
than large grains on top of small grains, this implies that 
#12 > #21, if type 1 grains are smaller than type 2 grains. 
Thus we arrive to the relations for the case of a granular 
mixture composed of small- rounded grains (type 1), and 
large-cubic grains (type 2), which gives rise to stratifica- 
tion: 



On the other hand, if the large grains are more rounded 
than the small grains, one expects that O22 < #11, and 
this type of mixture results only in segregation and not in 
stratification. The relation between the angles of repose 
is then 



721 < "22 < #11 < #12, 



(3) 



which is valid for small-cubic grains (type 1), and large 
rounded grains (type 2). 





FIG. 2. Different morphologies predicted by the present 
model, a, Stratification (S r — 0.2, ip r — 0.3, S u = 0.3, 
tp u = 0.3, U21 — 0.3, and r\2 = 0.1). Notice the "kink" formed 
by a pair of layers developed at the surface of the pile, sim- 
ilarly observed in Jl^,[L6|]. b, Strong segregation [5 r = —0.1, 
ipr = 0.3, 8 U — —0.3, npu = 0.6, u 2 i = 0.3 and ri2 = 0.1). 
Here the black grains are type 2, and grey grains are type 1. 
See Fig. [l] for the location of the morphologies in the phase 
space. 

Next, we generalize the model to the case of two types 
of grains 1 and 2 with different size, shapes or material 
properties. We deposit N\ and N2 grains on the first col- 
umn of the pile, and at a given time step, one rolling grain 
per species interacts with the sandpile surface, and the 
remaining rolling grains move downward to the nearest- 
neighbor column. We assume an overdamped situation 
where the rolling grains which do not interact with the 
surface achieve a constant convective velocity due to the 
collisions with other rolling grains. We define the gener- 
alized restitution coefficient and the generalized energy 
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barrier as r a p, and u a p with a,/? = 1,2, for the four 
different possible collisions, i.e. ri 2 is used in Eq. (]l|) 
if a rolling grain of type 1 collides with a static grains 
of type 2. Since the angle of repose of the pure species 
is a monotonic increasing function of the energy barrier 
u, and a monotonic decreasing function of the restitu- 
tion coefficient r, we can translate the relations @ for 
stratification and (|^) for segregation into relations for r a p 
and u a f3. Thus, we expect stratification for small-smooth 
grains (type 1), and large-rough grains (type 2) when 



"21 < Uu < "22 < "12 

ri2 < r 2 2 < rn < r 2 i, 



(4) 
(5) 



and, we expect segregation for a mixture of small-rough 
(type 1) grains, and large-smooth (type 2) grains when 



"21 < "22 < "11 < "12 

f*i2 < rn < r 22 < r 2 i. 



(6) 
(7) 
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■ "21 



- "12 

r\2 ~ 



~ "22 
f22, 



and 



K = "22 - "11 

S r = rn - r 22 . 



(8) 
(9) 



(10) 
(11) 



Further, we assume some values for S u and ip u (see Fig. 
|l|). We notice that ip u describes the difference in size 
of the grains, 5 U and 5 r are determined by the differ- 
ent shapes of the grains, and the material properties and 
asperities are described by tp r and S r . If grains 2 are 
rougher that grains 1, we have 8 r > 0, and if grains 2 are 
more smoother than grains 1, we have S r < 0. 

Figure [l] shows the resulting phase diagram and Figs. 
H and |H show the resulting morphologies. We find a re- 
gion of stratification when 8 r > (Fig. Ha), and a region 
of strong segregation for S r < (Fig. gb), as found in 
the experiments performed in Pq,E4,E5|. 
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FIG. 3. Morphologies predicted by the present model when 
the grains are only slightly different, a, When 5 r = and 
Vv < 0.3 (u2i = 0.1, ri2 = 0.4) we observe weak stratification 
with thin and irregular layers as observed in for small flux 
rate iVi = A*2=2, and b, mixing in almost all the pile plus 
weak segregation and some stratification with grains type 2 
at the bottom of the pile for large flux rate Ni — A*2=16. c, 
When 5 r > and ip r = 0.1 (1121 = 0.3, r\2 = 0.1) we observe 
the total mixing of the species. See Fig. [l] for the location of 
the morphologies in the phase space. 

We corroborate these predictions by investigating the 
different morphologies predicted by the model for the dif- 
ferent internal parameters. Since the general model has 
eight parameter, we reduce the number of them to be 
able to investigate the resulting phase diagram. We as- 
sume some relations between the parameters and define 
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FIG. 4. Wavelength of the layers (measured in pixel units) 
as a function of the flux of incoming grains (N\ + N2) for two 
different sets of parameters in the region of stratification. We 
find a linear dependence in agreement with conservation ar- 
guments , the coefficient of the linear relation depends on 
the internal parameters. The circles correspond to W21 = 0.1, 
ttu = 0.2, U22 = 0.6, U12 = 0.9, while the squares correspond 
to U21 = 0.5, uu = 0.9, «22 = 1.5, u\2 = 1.9. We take all 
the restitution coefficients equal to 0.2, and we use an equal 
volume mixture (iVi = iNfe). 

When some of the properties of the grains are very sim- 
ilar we find the additional morphologies shown in Fig. |^. 
When 6 r — and ip r < 0.3 we observe a weak strat- 
ification pattern (Fig. [5]a) but only for small flux rate 
Ni = N 2 = 2. The resulting layers are very thin and 
irregular, as can be seen in Fig. ||a. The negligible dif- 
ference in grain properties gives rise to identical angles 
of repose of the pure species, so that the kink, which is 
observed to give rise to the layers by stopping the rolling 
grains Jl6| , is very small. As a consequence, when we in- 
crease the flux of grains, the small kink is not able to stop 
the arriving rolling grains; the grains ride over the kink 
so that not segregation at the kink is observed. There- 
fore, for these particular parameters, the stratification 
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pattern dissapears upon increasing the flux. In Fig. ||b 
we show the results of our simulations where the same 
calculations of Fig. ||a are done, but for a larger flux of 
grains N\ = N 2 = 16. We see that at the 2/3 upper 
part of the pile the grains are mixed, and that at the 
1/3 lower part of the pile there is some reminiscence of 
stratification plus the segregation of the large grains at 
the bottom. In general we find that when Ni = N 2 ~ 8, 
the stratification disappears for these kind of parameters, 
a prediction that was recently confirmed by experiments 
pq| . Finally, when S r > and ip r = 0.1 wc observe the 
total mixing of the species for any value of the flux of 
grains (Fig. ^Jc). As in the previous case, the mixing 
might be due to the existence of a weak kink that is not 
able to stop and segregate the grains. 

We also study the wavelength A of the layers as a func- 
tion of the flux of grains (i.e., as a function of the total 
number N% + N 2 of deposited grains per unit time) . We 
find a linear increase of A (Fig. ^) as a function of the flux 
as is expected from a conservation argument |l^Jl^,^] : 

\<xNi+N 2 . (12) 

We find that the coefficient of the linear relation depends 
on the internal parameters of the model (Fig. ||). In gen- 
eral, we find that the larger the difference in energy bar- 
rier of the two species, the larger the wavelength of the 
layers, since the kink becomes steepest as the difference 
in angle of repose of the pure species increases. 
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